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Abstract 

The mean field Green function solution of the two-band singlet-hole 
Hubbard model for high-Tc superconductivity in cuprates (Plakida, 
N.M. et al., Phys. Rev. B51, 16599 (1995), JETP 97, 331 (2003)) 
involves expressions of higher order correlation functions describing 
respectively the singlet hopping and the superconducting pairing. Rig- 
orous derivation of their values is reported based on the finding that 
specific invariant classes of polynomial Green functions in terms of the 
Wannier overlap coefficients Vij exist. 

1 Introduction 

The two-band singlet-hole Hubbard model considered by Plakida et al. [T] 
for the description of the high-T^ superconductivity in cuprates in terms of 
Hubbard operators (HOs) provides the simplest consistent approach towards 
the incorporation of the essential features of these systems (strong antifer- 
romagnetic superexchange interaction inside the Cu02 planes, occurrence 
of two relatively isolated energy bands around the Fermi level, able to de- 
velop d^i-yi pairing) such as to describe simultaneously both the normal and 
the superconducting states within a frame which secures rigorous fulfilment 
of the basic principles of the quantum mechanics. The equation of motion 
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method for two-time Green functions |2I was successfully used to derive elec- 
tron spectra of the model P and to incorporate the superconducting state 
as well laH. 

The present paper reports rigorous results for the expressions of two 
higher order correlation functions which arise in the generalized mean field 
approximation (GMFA) solution of the Green function (GF) gj: {X^'^Xf) 
which describes the singlet hopping, and (X^'^Nj), which describes the ex- 
change superconducting pairing interaction. 

The idea proposed in ^ for the evaluation of the average (Xf^Nj) can be 
consistently generalized to yield power series expansions of both statistical 
averages. We find that the lowest order expressions of the two correlation 
functions, generically denoted henceforth (X^'^Qj) — where Qj is either 
or Nj, are obtained in terms of GMFA Green functions. 

This remarkable result follows from the mathematical properties of the 
Hubbard operators. These allow the definition of characteristic invariant 
classes of polynomial Green functions in terms of the Wannier overlap coef- 
ficients Uij, which are characterized by the property that, under the iteration 
of the equation of motion, the operator part remains invariant, while the 
polynomial degree in Uij is increased making thus possible consistent power 
series expansions in the small parameters Uij. 

2 Statement of the problem 

2.1. Hubbard operators. The Hubbard operators (HOs) X"^ = \ia){ip\ 
are defined for the four states of the model: |0) (vacuum), |cr) = If) and 
1^) = I i) (spin states inside the hole subband; in numerical calculations, 
a = ±1/2, a = — cr), and |2) = |||) (singlet state in the singlet subband). 

The multiplication rule holds X°'^X^^ = The HOs describing 

the creation/destruction of single states in a subband are Fermi-like ones 
and obey the anticommutation relations {X°'^ , Xj^} = 5jj((5^^Xf + (5,,aX7^), 
while the HOs describing the creation/destruction of singlets, spin and charge 
densities, particle numbers, are Bose-like ones and obey the commutation 
relations [X"^,XJ^] = 5ij{5f3^X'^^ — 5r,aX]^). At each lattice site i, the 
constraint of no double occupancy of any quantum state \ia) is rigorously 
preserved due to the completeness relation + X"" + Xf^ + XP = 1. 

The particle number operator at site i is given by 

X, =Xr + Xf + 2Xf . (1) 
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We define the Hubbard p-form of labels {a(3,'yri), 

rS''' = Y.^!,XfX]\ p=l,2,..., (2) 

where the meanings of HOs X"'^ and X'^^ depend on the context. The 
Wannier overlap coefficients Vij are small quantities rapidly decreasing with 
the intersite distance = |rj — rj| (see, e.g., [3 E] and references quoted 
therein). The nearest neighbour values i'i,i±a^/y = i/i ~ —0.14 and next 
nearest ones Vi^i^a^^ay = 1^2 — —0.02 considered in [T] are typical. 

2.2. Model Hamiltonian. The Hamiltonian of the model PP can be 
rewritten in terms of linear Hubbard forms describing hopping processes as 
follows 

1,(7 i 

+1^11 l^r^/ +IC22 2_^T^^i +/Ci2 2^2or(ri +ri/ ), (3) 

2,cr 2,cr 2,cr 

where the summation label i runs over the sites of an infinite two-dimensional 
square array the lattice constants of which, defined by the 

underlying single crystal structure. 

In Eq. El = Ed — II and E2 = 2Ei + Ueff, where id is the renormalized 
energy of a d-ho\e, /i is the chemical potential, while t/g// = A ^ Apd = Sp—Sd 
is an effective Coulomb energy corresponding to the difference between the 
hole energy levels for oxygen and copper. 

Keeping in mind that the lower label 1 refers to one-hole states, while the 
lower label 2 to singlet states, the quantities Kah = '^ipd.Kab are characteristic 
hopping energies for either inband (a = h) or interband (a 7^ 6; /C12 = /C21) 
transitions between the two bands of the model. Here tpd denotes the hopping 
p-d integral and Kab are numerical coefficients coming from hybridization 
effects between the holes and the singlets |lj. 

The translational invariance of the system gives 

which secures the hermiticity of the model Hamiltonian H. 

2.3. Mean field approximation. The quasi-particle spectrum and 
superconducting pairing within the model Hamiltonian (j2)) are obtained OE] 
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from the two-time 4x4 matrix Green function (GF) in Zubarev notation j2] 

- t') = {{XUt)\Xl{t'))) = -zOit - (5) 

where (■ ■ ■) denotes the statistical average over the Gibbs grand canonical 
ensemble. 

The GF (0) is defined for the four- component Nambu column operator 

K = {Xf Xf (6) 

The operator X]^ = {Xf" Xp X^) is the adjoint of Xj^. In (0), the 
superscript (T) denotes the transposition. Here and in what follows, due 
to translational invariance, the notation Qij points to the dependence of the 
quantity Q of interest on the distance rij = |rj — rj| between the position 
vectors of the lattice sites j and i respectively. 

The derivation of the GF within GMFA needs the knowledge of the fre- 
quency matrix, 

A,,= {{[X,,,HlXl}). (7) 

Direct calculations j3] show that the normal matrix elements of Aija 
contain one-site {i = j) GMFA hopping correlation function which result 
in 0{lCabVij) renormalizations of the energy parameters Ei and £"2, as well 
as two-site {i 7^ j) hopping generated higher order correlation functions 
bringing two distinct kinds of contributions to Aij^'- charge-spin correlations 
(which can be conveniently subdivided into charge- charge, (NiNj), and spin- 
spin, (SjSj) = {X^'^Xj"), correlations) and the singlet hopping correlation 
function (Xf^Xj^) (singlet destruction at site i followed by singlet creation 
at site j). 

There are three distinct matrix elements out of the eight normal matrix 
elements of Aija containing 

(a2,2a) - /Cn/^i,(Xf Xf ) (8) 

(Oa,aO) -/C22^i,(XfXf) (9) 

(a2,a0) -2a-/C2iz/i,(XfXf) (10) 

The only non-vanishing anomalous matrix elements of Aija are the 
hopping generated two-site contributions involving the higher order correla- 
tion function (Xf^Nj). This provides the exchange superconducting pairing 
mechanism originating in the interaction of an anomalous pair of particles at 
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a same site i but in different subbands = X^'^Xf^), with the surround- 
ing particle distribution at the neighbouring site j described by the particle 
number operator Nj, Eq. (|T)). The structure of the anomalous part of Aija 
is very special: 

(a2,a2) 2a-/C2iZ/ij(Xf iV,) (11) 

(Oa,o-0) 2(T-/C2iz/ij(XfiVj) (12) 

{a2,0cx) i(/Cii + /C22)/^.,(XfiV^) (13) 
(0a,a2) -l(/Cn + /C22)z/.,(XfiV,) (14) 
The other anomalous matrix elements are obtained by complex conjugation. 



3 Fundamental relationships 

From the spectral theorem |2], 

+ 00 

if duj 



^^^'^^^ = 2^ / 1^^ [((^f - {{XnQj)).~^e\ , (15) 



-oo 



where the labels ±ie, e = 0"*", refer to the retarded/advanced Green functions 
respectively. Since both and Qj are bosonic Hubbard operators, the 
thermodynamic factor in the denominator is 1 — e~^'^ and the two Green 
functions are defined in terms of the commutators of the two operators, i.e., 

{{xr{t)mm = -^e{t-mx^'{t),Q,{t')]) (le) 

for the retarded Green function and a similar definition for the advanced one. 

By differentiation with respect to t and use of Fourier transform, we 
get the following basic result for the two Green functions in the (r, u)- 
representation required by Eq. (fT^ (for the sake of simplicity, ±ie terms 
are omitted): 

{u; - E,){{X^'\Q,))^ =-/Cn ^(^Tl^:'-))- + E^^^ir'l^.))- + 

+/C2i$^2a((«f'^|g,)).- («ng,)).) (17) 
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Theorem 1 Let g2plJi = {{^2p-ii\Qj))^^ ^ generic notation of the exten- 
sions to Hubbard (2p— l)-forms ^ of the four Green functions which enter 
the r.h.s. of Eq. Then the recurrence relations hold, 



a2fla _ ^ij ^ij<T ^22^^ jj Pjjq ^ lC\i + /C22 cr2ficr 2/C11/C22 0cr,(T2 o\ 
92p-l — TP 1 \2^ i TP \2 92p+l I TP \2 92p+l U°J 



Qa, ct2 ^'^ija , ^^^^ij ^ij<T 2/C11/C22 a2,Qa , ^11 + ^22 0cr,<T2 /.qn 



^'^'-^ c^-2Ei (cu-2Ei)(cu-E2) (cu-2Ei)(c^-E2)^'^+' ^ ^ 

^2^-1 cu-(E2+A) [cu-(E2+A)](a;-E2) [cu- (E2+A)] (07-^2) ^'^+^ ^ ^ 

where the coefficients Mij^j and Pija, given in Table 1, are statistical averages 
following from equal time commutator terms. 

The proof is immediate if we write the equations of motion of the Green 
functions mentioned in the l.h.s. and iterate once. 



Table 1: Equal time commutators arising in the recurrence relations (jl8p - 
(nT\\ as coefficients of i^ff"^ and ufj 



2p-l 






2p 


Qi - 


g, = 


Ml,, 






P'- 








M" 

IJCF 




{xrxf) 


P". 


{xr{x^'-xf)) 


{xrxf) 


M'" 






pni 


-{xrxf) 









(x^xp) 


pIV 


{xr{xf-xf)) 


2(XfXf) 








pY 


{xrxf) 











pVI 


{xr{xf-xf)) 


2{XfXf) 



The p-form class invariance of the abovementioned Green functions is 
a straightforward consequence of the commutation and multiplication rules 
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satisfied by the Hubbard operators. It allows the derivation of the statistical 
averages (X^'^X'^^) and {X^'^Nj) as power series of the small parameters of 
the model ICuUij, JC22J^ij, and JC2ii^ij- 

While the coefficients of the odd powers of the series expansions are ob- 
tained in terms of GMFA Green functions, those of the even powers are still 
defined in terms of Green functions beyond GMFA. Prom the point of view 
of the practical implications, it is an academic question whether it would be 
possible to express them in terms of GMFA Green functions as well. The 
lowest order approximations provide the most important contributions to the 
observables. The question is whether the corresponding statistical averages 
quoted in the Table 1 are really significant or not. This point is discussed in 
the next section. 



4 Significant lowest order terms 

The results derived in the previous section show that the statistical averages 
are obtained as power series of Uij, with the contributions coming 
from the poles of the Green functions given by integrals of the form 

Imn{l^l,l^2) = Imni^l^^2) - ui + u > 2, m, n > 0, (22) 

i f div 1 1 

—oo 

The calculation of these integrals is standard: they are extended in corre- 
sponding complex half-planes z — {u!,'^z ^ 0), with half-circles at the three 
existing poles: z = 0, z = uji ± ie, z = 002 ^ is- The obtained contour 
integrals are calculated in two alternatives ways: using the residue theorem, 
and estimating them along the pieces of the involved contours. 

Retaining the lowest order contributions to (X^^Qj) only, we get: 

(XfQ,)=/2o(i?2)cf'°'^+/n(i?2.2i?i)cf'°'^+/n(£;2+A,£;2)cf''^'+0(4), 
cf'°-=-(2/Cni.,,-)M4..+(2/C22i/i,)M^„ 

cf''^'^=2<7(2/C2i^,,)^;, 
cf'''^2a{2K:2iiyij)Mg, 



ho{E2)-- 



IJC 



(5 El (1 - e-P^^Y 
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In{E2,2Eiy- 



PE2{E^A) A 1 - e-/3(^=2+A) A 1 - e-/^^^ • 



(24) 



Theorem 2 T/ie lowest order power series expansions of the correlation 
functions (X^'^Xj^) and (Xf^Nj) are obtained in terms of GMFA correla- 
tion functions as follows: 
For hole-doped systems: 

{X^Xf) ^ 2a^^(Xf Xf ); {XfN,) 2a^^(Xf Xf ). (25) 
For electron-doped systems: 

(Xf Xf ) ^ 2a^'^(Xf Xf ); (Xf X,) ^ 2a^'^(Xf Xf ). (26) 



Indeed, for hole-doped systems, the Fermi level lays in the upper singlet 
subband, such that we get the energy parameter estimates E'l ^ £'2 — —A. 
This yields 120(^2) ^ hi{E2,2Ei) ~ 0, while Jn(E2 + A,E2) - (2A)-i, 
therefrom Eqs. follow. 

For electron-doped systems, the Fermi level lays in the upper hole sub- 
band, and this yields the energy parameter estimates -Ei ~ 0, E2 — /S.. This 
results in 120(^2) ^ /ii(£2 + A, E2) ~ 0, while In (£2, 2Ei) ~ (2A)-\ hence 
Eqs. (I2S1) follow. 



5 Discussion of the results 

Corroboration of Eqs. ^ to (fT^ with the results stated in Theorem |21 show 
that, both in hole-doped and electron-doped cuprates, the dominant con- 
tributions to the singlet hopping and the exchange superconducting pairing 
are second order effects described by GMFA correlation functions. In (|25|1 
and there occurs a same small parameter, /C2it'ij/A, for the description 
of all the involved higher order correlation functions with, however, specific 
GMFA correlation functions. Thus, the singlet hopping proceeds hy i ^ j 
jumps of a particle from the upper energy subband to the lower energy sub- 
band. The anomalous superconducting pairing involves two spin states at 
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neighbouring lattice sites i and j, both with energies in that subband which 
crosses the Fermi leveL However, while both processes are given by small 
Oi^ufj) quantities, their consequences are quite different. 

The non-vanishing singlet hopping brings a small correction to the energy 
terms entering the normal part of Aija- Therefore, the decoupling ansatz 
(Xf ) ^ (Xf )(Xf ) = used in jT] does not substantially modify the 
general picture obtained for the normal state. 

On the other side, the derivation of the correct GMFA contribution 
to (X^'^Nj) is essential for understanding the pairing mechanism emerging 
within the model. Under the assumption of uniform hopping, Kn = K22 = 
K21 = K, Eqs. dnD-dUD and result in the well-known AFM ex- 

change interaction energy of the t-J model, J = At^/A, for spins on nearest 
neighbour sites and an effective hopping parameter t = tpdKui. 

In fact, the occurrence of different hopping coefficients in Eqs. (jHj)- (fT^ 
(ref. m, reported the values K22 ^ -0.477, Ku ^ -0.887, K12 = K21 ~ 
0.834) points to the existence of three asymmetric processes depending on 
the initial and the final energy subbands connected by the higher order cor- 
relation functions (Xf^Xj^) or (X^'^Nj) respectively. 

Finally, the present analysis strengthens the discussion in |^ concerning 
the unreliability of the approach of reference |7j towards the derivation of 
a GMFA expression for (X^'^Nj) based on the use of the Roth decoupling 
procedure which uncouples the Hubbard operators at a same site, 

(Xf iV,) = {X^'^XfN,) = {c,,c.,N,) ^ {{cUt)h,it')N,it'))). 

Therefore, the consequences following from this decoupling, namely that spe- 
cial interstitial excitations ("cexons") should appear and play an important 
role in the occurrence of superconductivity in cuprates, are artifacts of the 
procedure without actual physical meaning. 

6 Conclusions 

We reported a method for evaluating two higher order correlation functions 
describing respectively the singlet hopping and the superconductivity pair- 
ing within the two-time Green function approach to the solution of the two- 
band singlet-hole Hubbard model Q considered by Plakida et al. [H E] for 
the description of the physical properties of the high-Tc superconductivity 
in cuprates. Both in hole-doped and electron-doped cuprates, the dominant 
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contributions to the two processes are found to be second order effects de- 
scribed by GMFA correlation functions. 

The derived resuhs are rigorously established. For the two discussed 
processes, they rest on the occurrence of specific invariant classes of Green 
functions with respect to the operators of the Hubbard 1-forms entering the 
equation of motion (jl7|) . allowing therefrom powers series expansions in terms 
of the small Wannier overlap coefficients Uy. 

The singularity coming at the Fermi level from the thermodynamic factor 
(1 — e~^^)~^ is canceled by a corresponding singularity coming from the pole 
of the Green function at the Fermi level. From the point of view of the 
properties of the functions of complex variables, at = there arises a pole 
of the second order yielding the finite second order contributions found in 
the previous section. 

We have to remark that the existence of a non-vanishing commutator 
H] is essential. Since [A^,, H] = [Xf^, H] = 0, the method described in 
this paper cannot be used to the GMFA evaluation of the charge-charge and 
spin-spin correlation functions. 
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